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INTRODUCTION 

Dendritic crystal growth has been a subject of interest 
to physicists, metallurgists as well as mathematicians. 
The most common example of such a growth is the 
well-known ice-crystal. From a physicist’s perspective, 
dendrites constitute a relatively simple but important 
problem of pattern formation in non-equilibrium growth 
[3-5]. In metallurgy, dendrites are common to crystal 
formation in the manufacture of alloy when the growth 
rate exceeds some critical value. The literature on the 
subject is vast and reviewed in [3-6], as well as in our 
paper [1]. Here, we report a summary of the contents of 
our recent paper [2]; the reader is referred to it for more 
details. 

APPROACH 

In the first of a sequence of papers on dendritic crystal 
growth for weak undercooling [1], we derived asymp- 
totic equations for weak nondimcnsional undercooling 
A (non-dimensionalized appropriately, through a com- 
bination of latent and specific heat) for a dendrite that 
was asymptotically a parabola in the far-held. A Peclet 
number, P was introduced in accordance to 

A = V7P< p erfc (y/P) (1) 

which is clearly small for small A. Based on the length 
scale a, associated with the far-held parabola, a velocity 
scale r = 2 D P /ci was identihed, where D is the 
diffusion constant, a and ajV are used to nondimen- 
sionalize all lengths and times. We determined that if 
the initial deviations from an Ivantsov state (parabolic 
dendrite with a corresponding temperature prohle) are 
limited to an 0{ 1 ) region near the tip, then the dynamic 
evolution of the dendrite for the nondimensional time 
/ < P" 1 involves the 0(1) tip region only; in that re- 
gion, the temperature is harmonic to the leading order, 
with appropriate boundary and far- held matching condi- 
tions. It is to be noted that the derivation does not assume 
that the deviation from the Ivantsov state is small; only 
that it does not extend all the way to the far- held 0(P ~ 1 ) 
region. This tip-region dynamics was recast in terms of 
the evolution of the conformal mapping function from 
an upper-half C plane (C = £ + n?) to the exterior of the 


dendrite in the c -plane, where - = x + iy (See Fig. 1). 
This function t) was shown to satisfy the following 
nonlinear integro-differential equation for real C (i*e« on 


£-axis): 
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In the above, the nondimcnsional surface energy param- 
eter & is given by 


__ Cp jo Tm_ mj 

2a LP 

where do is the capillary length, c v the specihc heat, Tm 
the melting temperature of a planar interface, and L the 
latent heat. Further, in the above, k is physically the non- 
dimensional curvature and 1 L + a f is a four- fold surface 
energy anisotropy correction. Here ti the angle between 
the normal to the interface (pointing towards the melt) 
and the y axis, while 0 O is some fixed value denoting a 
direction along which surface energy is a minimum. 

Through a linearization of the equation for c(£,/) 
about some generally arbitrarily time dependent state, 
we were able to determine expression for growth rate 
of an initially localized disturbance in terms of the base 
state, through a Fourier-analysis, when the disturbance 
is far from the tip. In the special case of a base state that 
is steady and is close to an Ivantsov state, the expressions 
for the growth rate were in accordance to prior results 
[7]. Interestingly enough, it was possible to obtain the 
same results by analytically continuing the equations 
(2)-(8) to the lower-half complex C plane and carrying 
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out an asymptotic analysis for the linearized equations 
near singularities of 

It is to be noted that while the lower-half C plane does 
not correspond to any part of the physical domain, sin- 
gularities of c(C 0 approaching the real axis from below 
correspond to interfacial distortions. In particular, we 
found that according to the linearized dynamics, surface 
energy prevents an initially localized disturbance from 
remaining localized beyond a certain time. Arbitrarily 
small initial interfacial distortions (noise), representable 
by some singularity distribution in I m C < °> signifi- 
cantly affect the interface later in time when singularities 
of the associated zero- surface-energy problem approach 
or cross I w £ = 0 , even though surface energy locally 
smooths out all singularities in the linearized dynamics. 
The extent to which the zero- surface-energy singularity 
dynamics relates to growth rate and dispersion of dis- 
turbances for small non- zero surface energy was also 
uncovered. Hence, zero- surface-energy singularity dy- 
namics have both qualitative and quantitative impact on 
the physical predictions mentioned above. 

The relation between complex singularity dynamics 
and the evolving physical features of a dendrite tran- 
scends the restriction posed by linearized dynamics since 
a singularity of the conformal map in f w £ < 0 can re- 
sult in large interfacial distortions when that singularity 
approaches Im £ = 0 . In particular, if wc consider an 
isolated singularity C,(M of :< in the lower-half plane so 
that 

~ E 0 ( 0 (c-or 3 ( ] °) 

near < = C-( 0 . then lhe singularity is very close to 
the real axis, we can expect a comer pointed towards the 
crystal, with included angle (1 - ft) This is locally 
rounded off over a length scale determined by singu- 
larity distance from the real axis. The larger the |£o| 
(singularity ‘strength’), the larger is the impact region 
on the interface, mg E 0 determines the orientation of 
this distortion relative to the y-axis. The physical ef- 
fect of an isolated complex singularity corresponding 
to 3 = 1 (pole) is to create parallel sided indentation 
with width k E 0 and depth that scales as - In J 7 .., where 

>h = Im C-- , 

It is to be noted that the geometrical features at the 

interface associated with ( 10 ), as discussed above, will 
remain intact for a period of time, even when the actual 
singularity C,(0 > s smoothed out or replaced by acluster 
of other singularities, provided there is some intermedi- 
ate range: B* « |C - C40 1 « 1 for some s 311(1 
some set of real £ for which the behavior ( 10 ) persists. 

Prior work for dendrites [1], as well as by others 
on the mathematically analagous Hele-Shaw problem, 


shows that the zero-surface-cnergy dynamics preserves 
the form of the singularity - i.e. ,3 in (10) remains 
invariant with time; only its position C.40 ^nd its strength 
E 0 {t) evolve (except for a pole where E 0 is invariant). 
When 3 < 0, the form (10) is not invariant. Generally 
for an initial singularity of that kind, 

~ Ao(t) + Eo(t) (C - W0)"* (“) 

for ( sufficiently close to C- ■ Such singularity on the real 
axis does not introduce discontinuity in slope, except in 
non-generic cases- Aq = 0 just when Im £,* = 0- 

In this exceptional case, the comer is directed towards 
the melt, in contrast to 3 > 0 case, when it is directed 
towards the crystal. 

All singularities, regardless of their type, were shown 
to continually approach the real axis with time, though 
for 3 > j, they do not impinge the real axis in finite 
time-indeed they slow down significantly as they come 
close to the real axis. 

A point where = 0, but :,; is otherwise analytic, is 
referred to as a zero. A zero on the real axis corresponds 
to a zero-angled cusp on the interface that protrudes into 
the melt. Prior work, discussed in \\], shows that a 
zero remains invariant with time, when surface energy 
is neglected, i.e. the form 

~ 0 (C ~ Co(0) H-) 

remains invariant. The evolution equation lor £o(0> 
however, is found to be different from that of a singular- 
ity £,(7). In particular Co (A) ma y or ma y no1 a PP roach lhc 
real axis. For some set of initial conditions, a zero does 
impact the real axis in finite time. The mathematical 
solution ceases to be physically meaningful beyond this 
cusp- formation time. 

The connection between the dynamics in the ex- 
tended domain Im C < 0 and the physical features 

of an evolving dendrite, as described above, is partic- 
ularly useful, since there is strong evidence that the 
zero-surface-cnergy dynamics in the extended domain 
is well-posed [Sec [8]-[9] for evidence for the mathe- 
matically similar Hele-Shaw problem], in contrast to the 
interfacial evolution itself. In the latter case, the do- 
main is restricted to Im C = "This well-posedncss 
at the zeroth order mathematically justifies a systematic 
perturbation procedure in the extended complex domain 
to study how small but nonzero surface energy (with or 
without anisotropy) alters the zero- surface-energy dy- 
namics. The viewpoint we followed in [l]-[2], follow- 
ing the Hele-Shaw analysis with isotropy [9H10], is 
that the interfacial dynamics comes as a byproduct of 
the dynamics in the extended domain. 
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A necessary drawback to the above mentioned proce- 
dure is that now one must specify initial conditions in the 
extended complex domain Im £ < 0, which obviously 
cannot be done in an experiment where only the initial 
interface shape, up to some non- zero error, can be con- 
trolled. Connection to observed statistical features of an 
experiment can be made only by studying the statistics of 
an ensemble of complex-plane initial conditions, allow- 
ing for every conceivable singularity distribution, and 
with each member of the ensemble consistent with the 
given initial shape to within experimental error. Clearly, 
many different singularity distribution can result in the 
same approximate interfacial shape. However, an essen- 
tial precursor to such a statistical study is the thorough 
description of the dynamics of all possible forms for 
singularities in Im £ < 0. Once this is clarified, one 
can proceed with the statistical study for an ensemble 
of initial conditions. That such an approach may be 
useful is already demonstrated in [2], where we obtain 
dendrite coarsening results based on an ensemble of par- 
ticular singularities. However, in general, the analytic 
continuation of :(£,0) into Im £ < 0, corresponding 
to a general analytic initial shape, can be expected to 
contain natural boundaries and perhaps other singulari- 
ties that arc not isolated. Further, even the class of all 
possible forms of isolated singularities is too broad to 
study; only a small subset of possible initial conditions 
contains the specific classes of isolated singularities and 
zeros, as in (10) and (12), arc considered. Nonetheless, 
such isolated-singularity distributions do correspond to 
a range of interfacial distortion, when they come close 
to the real axis. For that reason, we believe that the sta- 
tistical features of the interfacial dynamics within this 
limited class of initial conditions are not very different 
from what is observed in experiment- with the additional 
proviso that a two-dimensional theory is applicable, at 
least insofar as scaling predictions. 

However, even within the class of possible initial 
singularities studied, there arc basic mathematical is- 
sues concerning the asymptotic matching of inner and 
outer regions in the complex plane (as the surface-energy 
parameter goes to zero) that remain unresolved. In car- 
rying out matching in the neighborhood of a singularity 
that is preserved by the zero- surface-energy dynamics, it 
is observed that the matching is necessarily sectorial-the 
inner solution does not match to the outer solution in ev- 
ery direction in the complex plane; it can be matched in 
a certain sector only. This is not a surprising result, since 
the steady dendrite problem is known to have the same 
features. However, unlike the steady problem where 
there are well defined global Stokes lines even beyond 


the immediate vicinity of an inner-region that determine 
local sectors of matching (See [ 1 1 ] for instance.), no ba- 
sic mathematical principle exists for the time-evolving 
flow. Only local Stokes lines, corresponding to local 
similarity solutions of the partial differential equations 
in the inner region, can be identified. There, we invoke 
a matching principle based on one used in the Hele- 
Shaw context [91 . The only direct evidence that such a 
matching principle is sound is our prior finding, in [1], 
that there is consistency between results from a Fourier 
analysis in the real domain and a complex singular! ty ap- 
proach involving inner-outer matching for the linearized 
problem. 


Further mathematical difficulties arise with initial 
zeros of since the full investigation of the dynamics 
at different stages is hampered by lack, in many cases, 
of either analytical or numerical solutions to a set of 
complicated partial differential equations in the com- 
plex plane. It is to be noted that the mathematical theory 
of nonlinear higher order partial differential equation 
in the complex plane is quite undeveloped. Progress 
in this case has been made, as in [8], with additional 
ansatz on the dynamics at intermediate stage[s]. There 
is no direct evidence that these ansatz are correct by 
themselves, though the careful numerical calculations 
of the interfaces themselves, for a sequence of compu- 
tations for decreasing surface energy, indirectly confirm 
the basic features of the analytic theory, both for the 
associated isotropic Hele-Shaw problem [10], and also 
for anisotropic Hele-Shaw and dendrite problems. This 
work will be reported in an upcoming paper. 


Despite the qualifiers above and the fact that our 
method necessarily requires a lengthy investigation of 
complex dynamics involving many kinds of initial sin- 
gularities with corresponding inner equations depend- 
ing on their distance from the real axis as well as the 
relative ordering of anisotropy and surface energy, this 
technique is the only one known for the fully nonlin- 
ear, time-evolving dendrite in the small -surface-energy 
limit. This limit is precisely the most difficult to explore 
computationally , since resolving small capillary lengths 
necessarily strains the capacities of computers. Further, 
even for cases where B is not small, the small surface 
energy limit cannot be avoided at large distances from 
the dendrite tip, where the curvature of an essentially 
parabolic keeps decreasing. 
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ISSUES ADDRESSED 

In [2], we continue our study of complex singularities 
initiated in [1] by including small but nonzero surface 
energy (0 < B << 1) in the nonlinear dynamics in the 
extended complex domain, generally taking anisotropy 
into account. The purpose of this paper is to address, 
partly or wholly the following important issues: 

1 . How does a non- zero B alter singularities described 
in (10)? Do the the alterations and modifications 
to the singularity stay confined to a small cluster 
around £,(/) ? Is there an intermediate spatial scale 
over which the behavior ( 1 0) is relevant as 6 — > 0 ? 
If so, is there a limitation on the order of | E 0 ( 0) | and 
the time for which this is so such behavior persists? 
How does anisotropy in surface energy come into 
play? 

2. What are the time and spatial scales over which 
surface energy effects become important to the real 
axis dynamics for a singularity corresponding to 
j > I ? Recall that according to zero-surface- 
energy dynamics, such singularities do not impinge 
the real axis in finite time, though it continually 
approaches it. 

3. For singularities corresponding to 0 < (3 < 
which are known, in the absence of surface energy, 
to impact the real axis in finite time (leading to 
comers at the interface), what are the relevant space 
and time scales associated with a small nonzero 6, 
when fm C.-(0 -» 0 

4. What can be expected about the growth rate of 
interfacial distortions associated with approaching 
complex singularities discussed above in (l)-(3) ? 
How does surface energy dissipation of weak sin- 
gularities, determine a cut-off in the growth rate ? 
How does anistropy affect the result ? It is often 
stated in the literature that intcrfacial distortions 
that point towards the crystal appear to remain sta- 
tionary in the laboratory frame. Is there a limitation 
on the time scale over which this is true ? 

5. What is the effect of anisotropic surface energy on 
an initial zero? Is there a ‘daughter singularity’ 
Q(0 that emerges from an initial zero 0(0), as 
for the isotropic Hele-Shaw problem [8]-[10] ? If 
so, how does anisotropy alter the structure of the 
cluster of actual singularities of that are centered 
at C40- 

6. How does the impact of 0(0 on the real axis affect 
the interfacial features? As with isotropic Hele- 
Shaw problem, can one expect the daughter singu- 


larity impact time to indicate when an actual in- 
terface will veer off from the corresponding zero- 
surface energy solution ? 

7. How are interfacial cusps, associated with a zero 
0(0 impacting the real axis in finite time, prevented 
by small surface energy effects ? One scenario is 
that small surface energy becomes important only 
when the interface becomes close to a cusp, i.e. 
when curvature of the zero- surface-energy solu- 
tion becomes large. The second is that the in- 
terface never comes close to cusp- formation be- 
cause it necessarily veers off from the correspond- 
ing zero-surface-energy solution significantly be- 
fore any Co(0 can impact the real axis. In the con- 
text of complex singularity dynamics of this paper, 
the two scenarios are distinguished by the question: 
does a daughter singularity 0(f) necessarily impact 
the real axis before the corresponding zero Co (/) ? 

8. How does a given disturbance, that may be as- 
sociated with many different complex singularity 
distributions cause 0(1) localized deviation in in- 
terfacial slope from a smooth background state, 
evolve in time. Is there a rescaling under which 
the equations remain invariant in the small surface 
energy limit ? What does such an invariance tell us 
about the dynamics ? 

9. How do surface energy and anisotropy modify or 
confirm the coarsening scenario that we proposed 
in [1]? The selection effect of surface energy on 
an ensemble of assumed singularities of different 
strengths is examined in [2], resulting in a predic- 
tion for the coarsening rale of \y\ lp for an interme- 
diate range of distances where \y\ is large compared 
to some inverse power of B. There is no necessary 
contradiction with the well known |y| 1/? coarsen- 
ing result [12], since that is valid for all sufficiently 
large distances from the tip. 
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Figure 1: Conformal map from upper-half £ plane to 
the dendrite exterior 
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